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Spin dynamics is considered in ferromagnets consisting of electron and nuclear subsystems inter- 
acting with each other through hyperfine forces. In addition, the ferromagnetic sample is coupled 
with a resonance electric circuit. Under these conditions, spin relaxation from a strongly nonequi- 
librium initial state displays several peculiarities absent for the standard set-up in studying spin 
relaxation. The main feature of the nonlinear spin dynamics considered in this communication is 
the appearance of ultrafast coherent relaxation, with characteristic relaxation times several orders 
shorter than the transverse relaxation time T2. This type of coherent spin relaxation can be used 
for extracting additional information on the intrinsic properties of ferromagnetic materials and also 
can be employed for different technical applications. 
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I. INTRODUCTION 



Spin systems can exhibit rather nontrivial dynamics when the magnetic sample is prepared in a strongly noncqui- 
librium state and, in addition, is coupled with a resonance electric circuit [1]. Due to the resonator feedback field, a 
coherent motion of spins can develop resulting in their ultrafast relaxation. However, the feedback field can organize 
coherent relaxation only when some initial mechanism triggers the process. A simple case could be the application 
of an external pulse at the initial time. If this pulse is sufficiently strong, the spin dynamics could be described by 
the Bloch equations. A more difficult, but interesting, situation is when no external pulse starts the process, but the 
latter develops in a self-organized way due to local spin fluctuations caused by their interactions. In such a case, the 
Bloch equations are not appropriate [1] and one has to resort to microscopic models. 

A theory of coherent spin relaxation in a system of nuclear spins inside a paramagnetic matrix has been developed 
[2,3], being based on a microscopic Hamiltonian with dipole interactions between nuclei. In the present paper we 
generalize this theory to include ferromagnetic materials. As far as ferrimagnets are often described as ferromagnets 
with an effective magnetization, our approach is applicable to ferrimagnets as well. In this way, our aim is to suggest 
a general microscopic theory valid for a wide class of materials, including those having long-range magnetic order. 



II. THEORY 

We consider a magnet consisting of electronic and nuclear subsystems interacting with each other by hyperfinc 
forces. The system of electrons possesses long-range ferromagnetic order. The subsystem of nuclear spins is prepared 
in a strongly nonequilibrium state, which can be done by dynamic nuclear polarization techniques. The sample is 
inserted into a coil connected with a resonance electric circuit. The general Hamiltonian describing a wide class of 
magnetic materials can be taken in the form 



in which 



is the Hamiltonian of electron spins, 



H — H e + H n + Hi n t , (1) 
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is the nuclear spin Hamiltonian, and 

H int = A e S 4 • I, + \ e E A t S i^i ( 4 ) 

i i^j a/3 

is the term corresponding to hypcrfine interactions. Here Jij is an exchange interaction; /i e = g e /J-B, with g e being the 

electron gyromagnctic ratio and /is , the Bohr magneton; the nuclear dipole interactions C"^ = ji n (j) a f3 — ^ n ?j n ij^j l r % 

contain fi n = g n ^N, with g n being the nuclear gyromagnetic ratio and /xjv, nuclear magneton, and = |ry|, riy = 
Tij /m , Vij = Vi— Yj ; the hypcrfine interactions consist of a contact part with a constant A and of a dipole-dipole part 

with A"? = (J, e (J, n i^afi — ^ n ?j n ij S j l r %'i t nc indices i and j enumerate electrons or nuclei according to the context, and 
a, (3 — x,y,z; Si is an electron spin operator while lj is a nuclear spin operator. The total magnetic field B is the 
vector sum of an external magnetic field in the z direction and of a transverse field Hi = H a + H in the x direction, 
consisting of an effective field H a of a transverse magnetocrystalline anisotropy and of a resonator feedback field H. 
The latter satisfies the Kirchhoff equation 

^ + 2 73 H + u? J* H(r)dr = -4th, ^ , (5) 

in which r\ is a filling factor; ui, resonator natural frequency; 73 = u/2Q is the resonator ringing width; Q is a quality 
factor; and M x — y J2i(^e < Sf > +fi n < If >) is the transverse magnetization, where the angle brackets mean the 
statistical averaging. 
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Employing the Hciscnberg equations of motion, we derive the time evolution equations for the following averages, 
related to the electron and nuclear spins, 

x = TT 5^ < S i~ > i z = TT < S i > > 

N e 1 Ne i 

--] v -E< / r>, -^E< J f>. (6) 

l z 

where AT e and AT„ are the number of electrons and nuclei, respectively, and S~ and I~ are the ladder operators. As 
the transverse variables x and u are complex, we need also the equations of motion for cither x* and u* or |x| and 
In this way, we obtain seven evolution equations, three for the electron variables x, z and | cc | 2 , three for the 
nuclear variables u, s, and |w| 2 , and the feedback-field equation (5). Although this is a rather complicated system 
of nonlinear equations, it can be treated by using the scale separation approach [2,3]. Details of this approach have 
been thoroughly described in Refs. [3-5]. 

Note first that if one invokes the standard semiclassical decoupling of spin correlators, assuming the translational 
invariance of the average spins, then some of the terms in the evolution equations become zero because of the properties 
of the dipolar interactions. The translation invariance of averages is equivalent to neglecting inhomogeneous local 
spin fluctuations. However, the latter are crucially important for the correct description of spin dynamics [1-3]. The 
inhomogeneous spin fluctuations can be retained by treating them as random local fields. Thus we come to the 
stochastic semiclassical approximation [3-5]. Then, using the method of the Laplace transforms, we may express 
the feedback field from Eq. (5) through the derivatives of spin variables and employ this relation in the evolution 
equations for x, z, \x\ 2 , and u, s, \u\ 2 . The latter spin variables can be classified into fast and slow with respect to 
each other by comparing their time derivatives. We keep in mind the following usual inequalities: 

\ji/uj e \ < i , \j 2 /co E \ < i , [ri/wjvl < i , |r 2/ w| < i , (7) 

in which 71 and Ti are the transverse attenuations for the electron and nuclear spins, 72 and T 2 are the longitudinal 
attenuations for electron and nuclear spins, respectively, and 

ue = (p-eHo - As) /h , u N = (n n H - Am) /% (8) 

are the electron spin resonance frequency and the nuclear magnetic resonance frequency, respectively; m being an 
average magnetization in the electron system. Another reasonable assumption is that the external magnetic field H n 
is stronger than the magnetocrystalline anisotropy field H a and that, similarly to (7), the inhomogeneous broadening, 
caused by local spin fluctuations, is smaller than the corresponding frequencies, so that 

\a e /uj E \ < 1 , \a n /u N \ < 1 , \i*/lo e \ < 1 , \T*/w N \ < 1 , (9) 

where a e = fi e H a /h and a n = ^, n H a /h are the anisotropy parameters and 7* and T* are the inhomogeneous widths 
for electrons and nuclei, respectively. Also, because the nuclear magneton is three orders smaller than the Bohr 
magneton, we have 

\Hn/He\ « 1 , ri/71 « 1 , r 2 / 72 « 1 . (10) 

Finally, we consider the case of a high quality resonator, having a large quality factor, and we assume that the 
resonator natural frequency is tuned close to the frequency of nuclear magnetic resonance, so that \A^/ujn\ <§C 1 and 
73 /u> <C 1 , where A jy = w — u>n ■ 

With these inequalities, we can classify all variables into fast and slow with respect to each other. Following the 
general scheme [3-5], we solve the equations for fast variables while treating slow variables as quasi-integrals of 
motion. Then, the found solutions for fast variables are substituted into the equations for slow variables and the 
right-hand sides of the latter equations are averaged over the periods of fast oscillations and over the random local 
fields. Introducing also the change of variables 

w = \uf- al + T } + 62 s\ ^V2^W^ m; (n) 

LO N LUN 

we come to the equations describing the slow nuclear spin variables 

rj q fill) 

-=r 25W -r 1 (.s-C), — = -2T 2 (l+gs)w, (12) 
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where 



r 2 w 



1 + 



(13) 



is the parameter of effective coupling of nuclear spins with the resonator 

For the relaxation times T\ = T^ 1 and T 2 
t <C T\, equations (12) can be solved analytically giving 



2 one usually has the relation T 2 <C T\. Therefore, for the times 



— ^-tanh 



t-t 
to 



- , w = \ —2- ) sech z 
9 V5 T o, 



t-t 

TO 



where tq is the collective relaxation time and to is the delay time, respectively, 



T 2 



TO = 



to 



^(l + gsoY+g 2 wl ' 
with sq = s(0) and wq = w(0) defined by initial conditions 
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T 2 - t (1 + gs ) 



T 2 +T a (l+gso) 



(14) 



(15) 



III. CONCLUSION 



We analysed the obtained solutions for the parameters typical of such ferromagnetic materials as EuO, EuS, EuSc, 
Li :E Fe3_ 2 ;04, Mn-rSbi^j;, NiMnSb, NiMnSi, Co 2 MnSi, and Co in the fee and hep phases [6,7]. Since ferrimagnets can 
often be modeled as ferromagnets with an effective magnetization [8] , our analysis is applicable as well to ferrimagnetic 
materials, such as MnFe 2 04. For these materials, taking as initial conditions so = —I, wq = 0, where / is a nuclear 
spin, we obtain the relaxation time tq = T 2 /g/, where g = n 2 fi e /2fi n , and the delay time to = To ln(2wjv/10 7 ). This 
gives the coupling parameter g ~ 10 4 , which, with T 2 <~ 10 -4 s, yields the relaxation time t ~ 10 -8 s. Because 
lon ~ 10 9 s _1 , we have the delay time to <~ 5 x 10 -8 s. 

The described regime corresponds to the ultrafast coherent spin relaxation, when during the time to + T o the initial 
strongly noncquilibrium spin polarization so = —I changes to its equilibrium value s(t) = / at t ^> to. The time 
to + To can be four orders less than the standard spin-spin relaxation time T 2 . This ultrafast coherent spin relaxation 
can serve as an additional technique for studying the spin-spin correlations in ferromagnetic materials, complementing 
other known techniques, such as neutron diffraction, light scattering, and nuclear magnetic resonance. The ultrafast 
relaxation mechanism can also find application in the important problem of fast repolarization of solid-state targets 
used in scattering experiments [9,10], as well as for fast switching devices in electronics and computing. 
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